Accepted by ApJ, 2010 February 9 

Preprint typeset using I^Tf^X style emulateapj v. 11/10/09 



O 

<N 

O 



O 

u 

Oh- 

6 

h 



> 

00 

<N 
<N 

<N 
O 
O 



A SYSTEMATIC SEARCH FOR X-RAY CAVITIES IN THE HOT GAS OF GALAXY GROUPS 

Ruobing Dong, 1 Jesper Rasmussen, 2,3 and John S. Mulchaey 2 

Accepted by ApJ, 2010 February 9 

ABSTRACT 

We have performed a systematic search for X-ray cavities in the hot gas of 51 galaxy groups with 
Chandra archival data. The cavities are identified based on two methods: subtracting an elliptical 
/3-model fitted to the X-ray surface brightness, and performing unsharp masking. 13 groups in the 
sample (~ 25%) are identified as clearly containing cavities, with another 13 systems showing tentative 
evidence for such structures. We find tight correlations between the radial and tangential radii of 
the cavities, and between their size and projected distance from the group center, in quantitative 
agreement with the case for more massive clusters. This suggests that similar physical processes are 
responsible for cavity evolution and disruption in systems covering a large range in total mass. We see 
no clear association between the detection of cavities and the current 1.4 GHz radio luminosity of the 
central brightest group galaxy, but there is a clear tendency for systems with a cool core to be more 
likely to harbor detectable cavities. To test the efficiency of the adopted cavity detection procedures, 
we employ a set of mock images designed to mimic typical Chandra data of our sample, and find 
that the model-fitting approach is generally more reliable than unsharp masking for recovering cavity 
properties. Finally, we find that the detectability of cavities is strongly influenced by a few factors, 
particularly the signal-to-noise ratio of the data, and that the real fraction of X-ray groups with 
prominent cavities could be substantially larger than the 25-50% suggested by our analysis. 
Subject headings: galaxies: active — galaxies: clusters: general — X-rays: galaxies: clusters 



1. INTRODUCTION 

Feedback and heating from active galactic nuclei 
(AGN) is considered a p rime candidat e for solving the 
"cooling flow" pro blem (iFabianl Il994f ) in the hot gas 
of ga l axy clusters ([B irzan et al. 2004; Dunn & F abianl 
[20061 : iRaffertv et all 120061 : IMcNamara fc NuTsenl [20071) 
groups, and giant ell iptical galaxies (| Jones et alT 2002: 
Machace k et al.l 12006), although the details of this pro- 
cess are not yet fully understood. Recent obser- 
vations using Chandra and XMM-Newton have pro- 
duced a large increase in the detection of X-ray sur- 
face brightness depressions ("cavities" or "bubbles") in 
many of these systems, interpreted as buoyantly ris- 
ing bubbles created by AGN outbursts. Studies of 
such cavities in individu al clusters, including Hydra A 
(IMcNamara et al.l 12000ft. Pe rseus dFabian et all 12 000). 
A2052 (jBlanton et al.l I2001D. A2199 (iJohnstone et al. 



l2002f h and Centaurus (Sander s fc Fabianll2002[ ). indicate 
that the outburst energy required to inflat e these cavi- 
ties would be sufficient to balance cooling (Bir zan et al.l 
120041 : IRaffertv et aTl l2QQ6). explain the l ack of gas cool- 
ing below T ^ 2 keV in cluster cores (jPeterson et al.l 
12001b iKaastra et al.l l2004). and repr oduce the bright end 
of the galaxy luminosity function ([p enson et al. 2003; 
Bower et al.ti2QQfl : iCroton et al.l l2006: Siiacki & Springe!) 



2006). 



Detailed studies of larger samples of X-ray bright ellip- 
tical galaxies and their surrounding cavities has further 
elucidated the possible role of AGN feedback for the ther- 
mal and morphological properties of the hot gas in these 
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systems. Bes t et aTl ([2005, 2006) combined observed cav- 
ity properties and cent ral radio powers of the systems in 
the lBirzan et al.l (j2004f ) sample with the inferred fraction 
of radio-loud ellipticals in the Sloan Digital Sky Survey, 
to show that the time-averaged heating rate by radio 
AGN in massive ellipticals can generally balance the cool- 
ing rate of the hot gas surrounding these galaxies. Using 
a sa mple of X-ray lumi nous ellipticals with identified cav- 
ities, All en et "aTl (j2006[ ) estimated AGN jet powers from 
cavity properties and snowed that these correlate tightly 
with the anticipated Bondi accretion rate of the central 
super massive black hole. Evidence that radio AGN may 
be able to affect their gaseous surroundi ngs also in lower- 
mass e llipticals has been provided by iDiehl fc Statlerl 
(2008a), who demonstrated that the amount of asym- 
metry in the hot gas morphology of ellipticals correlates 
with the central radio AGN luminosity, even down to the 
lowest detectable radio powers in relatively X-ray faint 
galaxies. 

These results all point to a close connection between 
central AGN radio outbursts and the creation of X-ray 
cavities in the surrounding gas. In terms of establish- 
ing the incidence and nature of such cavities, most 
work so far has focused on studying cavity properties 
and AGN interactions with the i ntracluster medium 
(ICM) within galaxy clusters (jBirzan et al.l 120041: 
IDunn fcJ Fabian 2004; IDunn et al.l 2005; Dunn & F abianl 



2006; Raffe rjy_eL al l 120061: "[St ernberg et al.l l2007t 



Birz an et al.l 120081 : IRaffertv et al.l 120081 : IDiehl et al.l 
2008). Results suggest that roughly 2/3 of X-ray 
bright cool-core clust ers harbor detectable cavities 
(jDunn fc Fabianl l2QQ6h . and that these cavities ap- 
pear to obey tight scaling relations between their 
size and projected clustercentric distance, offering a 
means of testing their nature and that of the inflating 
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mechanism (Die hl et al.l 120081 ). Comparable work on 
giant elliptical galaxies that do not repre sent central 
brigh t est cluster galaxies is s o far limited (j Allen et al.l 
2006; iMcNamara fe Nulsenl I2007D . but indicates a 
significantly smaller d etectable cavity fraction of ~ 1/4 
(jMcNamara fc Nulsenl 12 0071 ). However, studies of AGN 
heating and X-ray cavities within large samples of 
galaxy groups have largely been absent. Although AGN 
outbursts in groups are assumed to be smaller in scale 
and less energetic than those in clusters, they may 
play a more prominent role in the evolution of the host 
structure due to the shallower gravitational potential 
of groups. Estimating the incidence and properties of 
X-ray cavities in groups is therefore an important step 
toward understanding the role played by AGN outbursts 
for the evolution of baryons on group scales. 

In many X-ray bright clusters with high-quality data, 
X-ray cavities are prominent and can be easily identified 
as X-ray surface brightness depressions using visual in- 
spection. This can be augmented by radio data revealing 
ongoing AGN activity within the central cluster galaxy, 
or the presence of radio lobes coincident with the cav- 
ities. In this process, a number of factors could affect 
the detect ability of X-ray cavities however, including the 
position, orientation, and angular extent of the cavities 
jgnBlin fc Heinzl 120021 : iDiehl et al.l 120081 : iBriiggen et "all 
120091 ), along with observational details such as the sensi- 
tivity of the data. These issues are even more important 
in the group regime, where the lower intrinsic X-ray lu- 
minosities can further impede cavity identification. To 
search for X-ray cavities in groups, it is therefore impor- 
tant to consider additional methods that can aid simple 
visual inspection in identifying these structures. 

In the present work, we select a sample of 51 galaxy 
groups from the Chandra archive, in order to systemat- 
ically search for and characterize the presence of X-ray 
cavities in a large sample of galaxy groups. In addition 
to visual inspection of the raw X-ray data, we facilitate 
the detection process by employing two methods to first 
characterize the large-scale group emission, viz. unsharp 
masking and modeling of the surface brightness distribu- 
tion. A good handful of the groups in our sample have 
been previously reported to harbor detectable cavities. 
Here we find that at least 13 of the groups, and poten- 
tially as much as half of our group sample (26/51), show 
evidence of X-ray cavities. We quantify the size, group- 
centric distance, and relationship with the central AGN 
radio luminosity for all these structures, in addition to 
offering some considerations as to the preferred method 
for identifying cavities in these relatively X-ray faint sys- 
tems. 

The structure of this paper is as follows. In Section [2] 
we present the group sample and outline the data 
reduction process. Our approach to searching for 
cavities and establishing their properties are described 
in detail in Section [3l Section |4] presents the results 
of our Chandra analysis. To understand the efficiency 
of the two methods employed for cavity detection, we 
generated a set of mock data sets, the analysis of which 
is described in Section This is followed by a discussion 
of our results in Section [6] and conclusions in Section [71 
The cosmological parameters assumed in this paper are 
H = 73 km s" 1 Mpc" 1 , Q m = 0.27, and ^ A 0.73. 
Unless otherwise stated, all uncertainties are quoted at 



the 68% confidence level. 



2. SAMPLE SELECTION AND DATA REDUCTION 

The groups in this st udy were selected from the Chan- 
dra archival samples oflSxm et al.l (j2QQ9h (43 groups) and 
iRasmussen fc Ponmanl (|2007| ) (an additional 9 groups). 
To suppress the potential impact of instrumental arti- 
facts on our results, only groups for which the central 
bright regions are covered by a single ACIS CCD were 
considered. This excluded one group (A160) from the 
sample. Our final sample is listed in Tableland includes 
51 groups, ranging in distance from ~ 20-550 Mpc. A 
few of the groups have been covered by multiple Chan- 
dra observations. In those cases, we generally used the 
longest observation unless this had a large offset between 
the group center and the CCD aimpoint. We stress that 
the sample does not contain individual giant elliptical 
galaxies, nor any galaxy clusters, and that the two stud- 
ies from which our sample is drawn were not themselves 
designed for cavity studies. Hence, while the selection 
criteria employed in those two studies favor fairly undis- 
turbed X-ray bright systems, the sample should not be 
inherently biased toward systems with prominent cavi- 
ties. 

Standard data reduction and calibration was per- 
formed to all data sets starting from level one event 
files. Chandra Interactive Analysis of Observations soft- 
ware (ciao) v.4.1.1 and Chandra Calibration Database 
(CALDB) 4.1.1 were used in this work. New level one 
event files were created with the "acis_process_events" 
task in CIAO, including charge transfer inefficiency cor- 
rection, time-dependent gain adjustment, and screen- 
ing for bad pixels using the bad pixel map provided 
by the pipeline. For observations taken in Very 
Faint mode, additional background screening was per- 
formed. Grade/status filters were applied (excluding 
ASCA grades 1, 5, and 7) along with Good Time In- 
tervals filters to produce level two event files. Times of 
high background were eliminated based on lightcurves 
extracted in regions away from the extended source cen- 
ter. Images were then produced in the 0.3-2 keV band, 
chosen to optimize the signal-to-noise ratio at the typical 
group temperatures of our sample, using spatial bins of 
1 or 2 pixels depending on source extent and data qual- 
ity. Exposure maps were produced assuming a monoen- 
ergetic distribution of source photons at the peak flux 
energy (usually around 1 keV). The source image was 
normalized by the exposure map, correcting for the ef- 
fect of strongly variable exposure near the detector edges. 
We adopt a threshold of 1.5% of the maximum value of 
the exposure map, setting all pixels with exposure below 
this value equal to zero. 

3. DATA ANALYSIS 

To search for small-scale X-ray structure and identify 
potential cavities, we employ two methods: Modeling 
the surface brightness distribution of the groups using a 
two-dimensional (elliptical) /3-model and performing an 
unsharp masking procedure. 

3.1. Elliptical ft -Model Fitting 
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Fig. 1. — Illustration of elliptical /3-model fitting results for 
NGC5044 (top row), NGC 3402 (middle row) and ESO 552-020 
(bottom). From left to right: Input source image, fitted model, 
and residual image. 

With this method, we first aim to characterize the 
large-sc ale group emission by means of a n elliptical j3- 
model (Cava liere fc Fusco-Femia"nol Il976| ). We use the 
Sherpa package in CIAO to fit this model along with a 
uniform background to the exposure-corrected 0.3-2 keV 
images of all groups. Bright point sources were identi- 
fied visually and masked out in all fitting, and only data 
from the central CCD were considered. Free parameters 
of the model are /3, r c , ellipticity, position angle, and the 
normalization, in addition to the local background level. 

For each group, the best-fit model was subtracted from 
the input image to produce a residual image. For a 
good fit, the residual image will be flat almost every- 
where (modulo Poisson fluctuations), with any remaining 
structure revealing departures from the model, such as 
cavities. The top panel of Figure [T] shows an example of 
this, displaying the results for NGC 5044, the group with 
the largest number of source counts and the most promi- 
nent cavities within our sample. In this case, the two 
cavities are already visible in the original image but be- 
come much more prominent in the residual image. The 
middle panel of Figure [T] shows the case of NGC 3402, 
in which cavities are not clearly detected. Here the fit- 
ted model describes the group emission very well, and 
residual "structure" can be largely ascribed to Poisson 
fluctuations. Finally, images with low signal-to-noise ra- 
tio tend to have no cavities detected, as illustrated in 
the bottom panel of Figure HJ In this case, the relatively 
low number of counts precludes any robust conclusion 
regarding the presence of cavities. 

For some groups the model fitting fails to converge 
on sensible parameters, returning, e.g., a best-fit core 
radius smaller than one pixel or a model eccentricity 
higher than 0.99. The relevant systems are 3C442A, 
A1139, A1238, A1692, A1992, A2092, A2462, NGC 507, 



RBS461, RXJ 1206-0744, and UGC842. We ascribe the 
failure of the fitting approach mainly to the fact that 
many of these data sets have relatively low signal-to-noise 
ratio, coupled with the possibility that some of those 
groups are simply not well described by an elliptical f3- 
model. We will return to these issues in Sectional 

3.2. Unsharp Masking 

Another way to smooth the large-scale group emission 
and test for structure on smaller scales is unsharp mask- 
ing. This approach has been successfully employed to 
uncover faint feat ures at fine spatial d etail in deep cluster 
X-ray data (e.g., Fabian et al.l 120061 ). Here we perform 
this procedure by smoothing the exposure-corrected im- 
ages using the "aconvolve" task in CIAO. For each group, 
the data are first smoothed using a wide Gaussian ker- 
nel that preserves the overall morphology of the emission 
but erases small-scale structure. A separate smoothing 
is performed on smaller scales to suppress (partly noise- 
induced) pixel-to-pixel variations while preserving struc- 
ture on the likely scale of any cavities in the data. The 
latter image is then divided by the former, with the re- 
sulting quotient image acting as an analog of the residual 
image from the model fitting approach. In principle, op- 
timal choices for the characteristic smoothing scales vary 
among the groups. In practice, however, we found that 
large and small scales of 10-30 and 2-5 pixels, respec- 
tively, generally produced the most visually compelling 
results, with the smaller values generally preferred for 
more distant and compact systems. Figure [2] illustrates 
this method applied to NGC 5044. The Figure shows the 
original image, the image smoothed by narrow and wide 
Gaussians of a = 3 and 30 pixels, respectively, and the 
result of dividing the former by the latter. The result can 
be directly compared to the residual image from model 
fitting shown in the top row in Figure [TJ 

To illustrate the impact of different choices for the 
small smoothing scale, Figure [3] shows the quotient im- 
age for NGC 5044 with six different small scales employed 
along with a fixed large-scale kernel of a = 30 pixels. 
All the quotient images clearly reveal X-ray depressions 
at the positions of the cavities, but it is easier to vi- 
sually recognize the cavities as such in the images pro- 
duced with relatively narrow small-scale kernels: The 
cavities in these images are more obvious and have a 
higher contrast with their surroundings. On the other 
hand, when the small scale becomes larger than 8 pixels 
(i.e. > 1/4 of the larger scale), distinguishing the cavities 
from other brightness depressions becomes more difficult 
as the sharp contrast between the cavities and their im- 
mediate surroundings is lost. However, we emphasize 
that NGC 5044 is by far the highest-flux system in our 
sample, and that the choice of small smoothing scale is 
less straightforward for significantly fainter systems. As 
such, a major drawback of this method compared with 
surface brightness modeling is the need to make appro- 
priate choices for the smoothing scales. 

3.3. Cavity Identification and Characterization 

Based on the images described in Sections I37T1 and \3?2\ 
visual inspection for X-ray cavities was performed for all 
groups. Cavities were mainly identified in the residual 
images from model fitting, but the appearance of the un- 
smoothed exposure-corrected images was also taken into 
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Fig. 2. — Illustration of unsharp masking for NGC 5044. From left to right: exposure-corrected image, small-scale smoothed image 
(Gaussian of a = 3 pixels), large-scale smoothed image (<r = 30 pixels), and quotient image. 

forms a ring around the group center. This is likely an 
artifact of the fitting procedure, owing to the presence of 
a strong central and extended excess above the best-fit 
/3-model. 

We emphasize that cavity detection was done solely on 
the basis of X-ray data, and that the presence or mor- 
phology of central radio emission in a group was not con- 
sidered for this purpose. Furthermore, as indicated by 
comparison with other existing studies (see Section HJ, 
we have likely been conservative in our classification. It 
is, therefore, conceivable that groups in the P-sample 
may generally have a high probability of containing cav- 
ities. 

In addition to identifying cavities we also provide quan- 
titative estimates of their properties. In each case, the 
cavity center was first defined as the centroid of the X- 
ray surface brightness depression. In a few cases where 
the deepest depression point of the cavity in the resid- 
ual image from model fitting differed significantly from 
the apparent geometrical center, we considered the cav- 
ity center to be the midpoint between the two locations. 
The location of cavity boundaries were identified as sharp 
drops in X-ray surface brightness, equivalent to where the 
pixel counts in the residual images become negative. On 
this basis, the cavities can be generally viewed as ellipses, 
with a major axis a in the tangential direction (perpen- 
dicular to the line joining the cavity with the group cen- 
ter) and a minor axis b in the radial direction. The cavity 
distance D was defined as the distance between the cav- 
ity center and the group center, with the latter defined 
as the location of the central peak of the extended X-ray 
emission. Since the accuracy of these somewhat subjec- 
tive measurements generally depends in an unquantifi- 
able manner on factors such as the signal-to-noise ratio 
of the data, we do not provide uncertainties on these re- 
sults, and we stress that they should all be considered 
approximate. We do note, however, that the inferred 
properties of the cavities in the P-sample are generally 
more uncertain than those of the C- sample. 

4. RESULTS 

Based on the above criteria, 26 out of the 51 groups in 
our sample were identified as harboring certain or possi- 
ble X-ray cavities. These systems are evenly split among 
the C- and P-samples. In Figure |4] we show the input 
image, the residual image from model fitting, and the 
quotient image from unsharp masking for all groups in 
the C-sample. Figure [5] shows the corresponding results 
for the P-sample. For further comparison, we also show 



Fig. 3. — The effect of choosing different widths of the small-scale 
Gaussian kernel in unsharp masking for NGC 5044. From left to 
right, top to bottom: quotient images with a = 2, 4, 6, 8, 10, and 
15 pixels (the large-scale kernel has a = 30 pixels). 



consideration. In most cases the impression gained from 
these two images agree. In the cases for which the fitting 
approach did not converge on reasonable model values 
(typically for sources of low signal-to-noise ratio), cavities 
were identified visually based on the raw images, aided 
by the quotient images from unsharp masking. For the 
groups judged to show evidence of cavities, differences 
in signal-to-noise ratio and in the contrast between the 
cavities and their surroundings led us to further subdi- 
vide the sample, leading to the following classification 
scheme: 

Groups with certain cavities (denoted the "C" -sample 
in the following): Cavities are clearly detected upon vi- 
sual inspection of the residual images, and their presence 
is at least indicative in the raw and unsharp masked im- 
ages as well. These cavities always present a high con- 
trast with their surroundings, including the presence of 
a bright rim. 

Groups with possible cavities ("P" -sample), fulfilling 
one of two conditions: (1) Cavities are apparent in only 
one of the residual and raw images. (2) Both the residual 
image and the raw image show a hint, if not conclusive, of 
the presence of cavities. Cavities in this category always 
have a low contrast with their surroundings. 

Groups with no cavities ("N" -sample): Neither raw 
nor residual images show any visually obvious evidence 
of cavities. This category also includes cases in which 
the surface brightness depression in the residual image 
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Fig. 4. — The 13 groups in the C-sample with clearly identified cavities. For each group, left panel shows the exposure-corrected image, 
middle panel the residual image from model fitting smoothed by a Gaussian of a = 2 or 3 pixels, and right panel the quotient image from 
unsharp masking. Ellipses overlaid on the residual images outline identified cavities. Horizontal bars mark a physical scale of 10 kpc. 
Where appropriate, vertical bars are used to indicate a different scale. 
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Fig. 5. — As Figure 2] but for the 13 groups in the P-sample with tentative evidence of cavities. 
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Fig. 6. — As Figure [4] but for the 25 groups in the N-sample without clearly identifiable cavities. The residual images in this figure have 
not been smoothed. 



the images for the N- sample in Figure [H 

Among the 26 groups in the combined C- and P- 
samples, about half (8 in the C- sample and 6 in the 
P- sample) show evidence of containing two cavities, 
and these are always at symmetric positions relative 
to the group center. We did not identify any groups 
showing evidence for more than two cavities as oth- 
erwise reported fo r a number of galaxy clusters (e.g., 
iSanders L Fabianl 120071 : IWise et al.ll2007[ ). At least nine 
groups in our sample have been ide ntified in previ- 
ous work as hosting cavities: IC 1262 (jTrinchieri et al.l 



2007h . A262 (Bir zan et al.ll200|), HCG62 dMorita eTal] 
2006), NGG507 ( Kraft et~aTl l200i IAllend~al.l 12006ft 
NGC 5846 (| Allen et all 120061). N GC741 (jJetha et al.l 
l2007h. NGG4325 (Russell etS] l2007h . NGC 5044 



(iGastaldello et aLll2QQ9h . and NGC5Q98 (jRandall et al.l 
2009). Eight of these were classified as belonging to the 
C- or P-sample in our analysis. The only exception is 
NGC 4325, which shows weak hints of the presence of 
cavities but without satisfying our criteria for inclusion 
in the P- sample. The i dentifi cation of cavities in this 
system bv iRussell et al.l (j2007l ) was based on a detailed 
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Bayesian model-fitting approach which clearly exceeds 
the level of sophistication employed in our analysis. Fur- 
thermore, two of the above nine groups (NGC 507 and 
NGC 741) were here classified as belonging to the P- 
r at her than C-sample. These comparisons confirm our 
suspicion that we have been fairly conservative in our 
classification. This might be anticipated given that we 
are not incorporating, e.g., radio data to aid in cavity 
identification. 

We list the /3-model fit results and cavity measure- 
ments (where applicable) for all the groups in Table [2j 
along with the total number of exposure-corrected 0.3- 
2 keV photons from diffuse emission within the model 
fitting region considered for each group (typically the 
central ~ 5' x 5 ; ). In agreement wit h many previous 
studies (e.g. JOsmond fc Ponm an 2004), we note that the 
mean f3- value derived for the sample is f3 = 0.47, con- 
siderably lower than the typical value of ^ 2/3 se en in 
clusters (jArnaud fc Evrardlll999l : iMohr et al1 fl999). We 
also tested for a relationship between j3 and core radius 
r c for all three subsamples of groups, as illustrated in 
Figure [71 No clear correlation between these two param- 
eters is evident, neither for the individual subsamples nor 
for all groups combined. In addition, mean values of the 
two parameters do not show any clear variation among 
the subsamples, being (r c = 7.9 kpc, f3 = 0.45) for the 
C-sample, (r c = 8.8 kpc, f3 = 0.46) for the P-sample, 
and (r c = 6.5 kpc, f3 = 0.49) for the N-sample. 

From the results in Table [2j we can explore any cor- 
relations among the sizes and groupcentric distances of 
the identified cavities. Figure [8^a) shows the relationship 
between the cavity size in the tangential direction a and 
radial direction b for the C- and P-samples. Remarkably, 
these two quantities appear strongly coupled, despite the 
potentially considerable systematic uncertainties associ- 
ated with their estimation. Motivated by the appearance 
of the data in Figure [8j we fitted linear relations to the 
data in log-log space to quantify the observed cavity be- 
havior. For this, we used t he orthogonal regression ap- 
proach of llsobe et al.l (| 19901 ), since there are (unquantifi- 
able) errors on both parameters, but we note that results 
are consistent with those obtained using various ordinary 



least squares methods (e.g., Y vs. X, X vs. Y, bisector). 
For the data in Figure Ufa), we then find 

loga = (1.00 ± 0.03) log6± (0.22 ±0.02), (1) 

for a and b in kpc and with a high correlation coeffi- 
cient of 0.95. This shows that cavities maintain broadly 
similar shapes, with a mean ratio of a/b ~ 1.7 that is 
roughly constant for all cavity sizes. Cavity sizes are 
also strongly correlated with their projected distance D 
from the group center. Figure [8fb) and (c) show a and 
b as functions of D (also in kpc), with best-fit linear re- 
lationships of 

loga = (0.86 ± 0.04) \ogD ± (0.20 ± 0.03), (2) 

for a correlation coefficient r = 0.94, and 

log 6= (0.87 ±0.03) log D- (0.03 ±0.02), (3) 

for a coefficient r = 0.88. As illustrated by the dotted 
line in Figure [5Jb), equation (j2j) implies that cavities are 
enclosed by a cone of roughly constant opening angle 
~ 60° as they rise and expand within the surrounding 
medium. 

In Figure[8j the fits of equations (P)-© are represented 
by solid lines. For comparison, we also show by dashed 
lines our corresponding fits to the results for the cluster 
sample of Birzan et al. (2004; their table 3), namely 

loga = (1.07 ±0.07) log b ±(0.10 ± 0.07), 
loga= (0.92 ± 0.07) log L> ±(0.13 ± 0.07), 
log b = (0.85 ± 0.06) log D ±(0.04 ± 0.06). 

In addition, Figure [8fb) furt her shows the relation be- 
tween a and D found by iDiehl et all ([20081 ) for an even 
larger cluster sample. Despite the somewhat subjective 
way cavity sizes are estimated by different authors, there 
is generally remarkable agreement between these three 
studies and hence between results for our groups and 
more massive clusters. In particular, the slopes of equa- 
tions H])-(j3]), along with two out of three intercepts, are 
statistically consis tent at the 1-cr level with the results of 
iBirzan et all (|2004fl . 

We next tested for a correlation between the 1.4 GHz 
radio luminosity Li.4gh z of the central brightest group 
galaxy (i.e., that of any central radio source associated 
with the galaxy itself rather than that of any radio lobes 
coinciding with the X-ray cavities) and the detection 
of cavities. Values of Li4GHz were extracted from the 
NRAO VLA Sky Survey (NVSS), and are listed in Ta- 
ble [2j If data were missing at 1.4 GHz but available 
at a nearby frequency, we extrapolated to 1.4 GHz as- 
suming a power-law spectrum of index unity, S oc v~ x . 
For groups with NVSS coverage but with no radio flux 
listed in NED, we used the 99% completeness flux limit 
in the NVSS survey (3.4 mJy; ICondon et all I1998D to 
estimate an upper limit to Li.4GHz- Finally, if a source 
was not covered in NVSS or if the 1.4 GHz emission could 
not be unambiguously asso ciated with the cent ral group 
galaxy in NVSS or FIRST (jBecker et al.lll995l ) data, we 
left I/1.4 GHz as undetermined. 

Figure [9] shows the resulting distribution of Li.4gh z 
for the different group subsamples. For groups with up- 
per limits to Li.4GHz, we have assumed this upper limit 
when displaying these results, but we note that the re- 
sults would not be substantially affected by excluding 
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Fig. 8. — (a) Cavity major axis a as a function of minor axis 
b for the C- and P-samples. Dotted line represents equality, (b) 
Major axis as a function of projected cavity groupcentric distance 
D. Dotted line illustrates the (arbitrarily normalized) expectation 
for cavities enclosed by a cone of constan t opening angle fl , and 
dot-dashed line shows the fitted relation of Diehl et al. (2008|). (c) 
Minor axis as a function of D. Solid lines in all plots show the best- 
fit linear relations for our group sample, equations while 
dashe d lines represent our fits to the cluster results of lBirzan et al.l 
(12001 ). 

those groups. There is no strong link between the pres- 
ence of detectable cavities and the current radio power 
of the central group galaxy. For example, the geometric 
mean of L\ 4 GHz is very similar for the three subsamples, 
10 22 - 44 W Hz" 1 (C-sample), 10 22 58 W Hz" 1 (P), and 
10 22 - 49 W Hz -1 (N). Again, excluding groups with upper 
limits to I/1.4 GHz does not change this conclusion. We 
also note that no strong trend is seen between projected 
cavity distance D and Li.4GHz, with the data showing a 
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Fig. 9. — Distribution of central 1.4 GHz radio luminosities for 
the different subsamples. 

small correlation coefficient of 0.4. 

Finally, to search for evidence of a 1.4 GHz radio 
plasma at the location of the cavities, we have compared 
the X-ray and radio morphology of the groups in the C- 
and P-samples, using radio data extracted from NVSS 
and FIRST. Since the cavities are generally present on 
spatial scales well below the resolution of NVSS data, 
comparison between Chandra and NVSS data generally 
provides no clear evidence of any morphological similar- 
ities. However, 12 of our groups have higher-resolution 
FIRST data in NED, and of those 12 systems, we show 
the six groups belonging to the C- or P-samples in Fig- 
ure |T0j No clear association between the X-ray cavities 
and any radio lobes is seen, however. It remains a pos- 
sibility that the cavities coincide with radio emission at 
frequ encies well below 1.4 GHz (e.g., iGiacintucci et al.l 

HqqI). 

5. TESTING CAVITY DETECTION WITH MOCK 

DATA 

To test the performance of the adopted methods for 
cavity detection and better understand the biases inher- 
ent in either method, we next generated a set of mock im- 
ages designed to match typical Chandra exposures within 
our sample. These images were then analyzed in a man- 
ner similar to that of the real data. This Section de- 
scribes the generation and analysis of these mock data. 

5.1. Model Setup 

The mock images were generated by first modeling the 
group emission as an elliptical /3-model. Cavities were 
superposed as circular depressions in local surface bright- 
ness, and a spatially uniform background was added. 
The following sets of model parameters were considered: 

• $ = {0.35, 0.50, 0.65}. 

• Core radius r c = 10 pixels. 

• Total number of source counts N tot = {5,000, 
15,000, 30,000}. 

• Cavity "strength" (the local surface brightness de- 
pression factor): {1.5, 2.0, 2.5}. 
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Fig. 10. — Comparison of X-ray and radio morphology for the 6 groups in the C- and P-samples with available FIRST radio data. For 
each group, left panel shows the exposure-corrected image, center panel the residual image, and right panel the FIRST 1.4 GHz image. 
Contours extracted from the radio data have been overlayed on the residual images. The first five groups belong to the C-sample, while 
UGC 5088 belongs to the P-sample. 



• Projected cavity distance D from X-ray peak in 
pixels: {15, 35}. 

This choice of parameters was generally motivated by 
the results for our real sample. The adopted j3- values 
span the range of fitted values covered by all but four 
of our groups. To keep the number of output images at 
reasonable levels, r c was fixed for all mock data. This 
should not introduce any substantial systematic bias in 
the results, since Figure [7] shows that r c does not vary 
systematically with f3 for our groups. Values of iVtot 
were chosen to range from the median of the groups with 
no detected cavities (~ 5000) to a value below that of 
HCG 62, one of the groups with the most prominent and 
obvious cavities. For the cavity radii, we assumed R = 
0.6 D, in rough agreement with equation ((2]). 

A total of 54 mock images were generated to span all 
combinations of the above parameters. For each image, 
a uniform background level was added, drawn from a 
Gaussian distribution with a mean equal to the approx- 
imate 0.3-2 keV value for a typical 50-ks ACIS-S expo- 
sure (« 0.1 counts arcsec -2 ), and with a = 0.2 x mean. 
The minor-to-major axis ratio of the /3-mode l was drawn 
from a Gaussian with mean 0.8 and a = 0.12 (Moh r et al.l 
1995), but was restricted to > 0.65. Position angles of 
the /3-model and of the cavities themselves were indepen- 
dently drawn randomly from a uniform distribution. The 
resulting surface brightness model was then convolved 
with the Chandra point spread function, derived using 
the "mkpsf" tool in CIAO at a photon energy of 1 keV 
and at the detector position of the ACIS-S3 aimpoint. 
Finally, Poisson noise was added to the mock images. 
Other instrumental effects (or point sources) were not in- 
cluded, since the model fitting to the real data was per- 
formed on exposure-corrected and point-source excised 
images. 

Admittedly, the model setup is a rather simplistic one; 
we assume the number of cavities in each system is al- 



ways two, they represent similar (local) brightness de- 
pressions, appear circular, and are at the same distances 
from the center of the X-ray emission. As such, orienta- 
tion effects were not taken into account (e.g., the cavities 
were assumed to be in the plane of the sky). However, 
since the purpose of this exercise was simply to test and 
understand our ability to recover cavities in the real data 
under reasonably favorable conditions, a detailed explo- 
ration of the full cavity parameter space would be beyond 
the scope of this work. 

5.2. Cavity Detectability 

We applied both the methods described in Section [3] 
to the mock images in a test of the veracity of either 
approach. Figure [IT] shows two example mock images, 
with A^ot =30,000 and 5,000, respectively, along with 
the results of our two methods. For the brighter source, 
the input cavities are easily recognizable in both output 
images. For the fainter case, however, there is some am- 
biguity as to the apparent location of cavities in the two 
output images. We will return to this issue below. Fig- 
ure [12] shows the distributions of best-fit values of core 
radius and j3 obtained from surface brightness modeling 
of all the mock images. Encouragingly, all distributions 
are peaked at the input values, lending credibility to the 
use of this method for describing group emission on large 
scales. For example, even though the mock sources are, 
on average, fainter than our real ones, there is no over- 
lap between the resulting distributions for the three input 
values of j3. Figure [13] summarizes the results of searching 
for cavities in both the output mock images, the quotient 
images, and the residual images. While the rate of "cer- 
tain" cavity recovery is comparable for the unprocessed 
and quotient images, using the residual images is clearly 
more efficient, with the success rate increasing by a factor 
of two (from 14-17 to 30 out of the 54 data sets). 

We found that all four controllable input parameters 
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Fig. 11. — Two examples of raw and processed mock images: 
Top panel shows a model with total source counts Ntot =30,000, 
(3 = 0.65, cavity strength 2.5 (prominent cavities) and cavity dis- 
tance D = 35 pixels. Bottom panel: TV-tot = 5,000, (3 = 0.65, 
cavity strength 1.5 (weak cavities) and D = 35 pixels. From left to 
right: Input source model, raw mock image, residual image from 
the /3-fitting method, and the quotient image from unsharp mask- 
ing (both slightly smoothed). 
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Fig. 12. — Distribution of fit results for (top) core radius r c and 
(bottom) p for the mock images. Vertical dashed lines represent 
the mock input values. 



(/?, iVtot? cavity strength and distance) influence the iden- 
tification of cavities, even when applying the relatively 
successful model fitting approach. This is illustrated in 
Figure HU which outlines the impact of changing the var- 
ious model input paramet ers, again using the classifica- 
tion described in Section 13.31 In particular, the value 
of TVtot emerges as a crucial factor for cavity recovery. 
Among the 18 mock images with iV tot = 5,000, only eight 
were identified as hosting possible cavities in the resid- 
ual images, while for N tot =30,000, the corresponding 
value is 17 out of 18. Remaining input parameters also 
play a role. For example, at fixed iVtot, cavity recovery is 




No Cavities Detected 



Possible Cavities 



Certain Cavities 



Fig. 13. — Number of mock data sets identified as harboring 
certain, possible, or no cavities based on the unprocessed mock 
images (empty bars), quotient images (grey), and on the residual 
images from model fitting (black). 

more efficient for j3 > 0.5 compared to j3 = 0.35, and for 
higher cavity strength and smaller distance. Another les- 
son from the mock data was that Poisson noise can easily 
mask one of the two cavities in X-ray faint sources. In 
other words, visual identification of only one cavity in 
unprocessed data does not necessarily exclude the pres- 
ence of a second cavity. Finally, we note that we found 
no clear correlation between the identification or input 
properties of cavities and various statistical moments of 
the photon number distribution in the residual images. 

Naturally, we cannot with confidence exclude the pos- 
sibility that a few cavities have been falsely identified in 
the Chandra data. This problem necessarily also plagues 
all other similar studies to some extent, but could be 
particularly acute for the present work since we have not 
incorporated radio data or other means to aid in cav- 
ity detection. However, we have reasons to believe that 
this issue is not important within our present analysis. 
First, as already discussed in Section HJ our comparison 
to existing results suggests that we have been rather con- 
servative in our cavity identification. Second, none of the 
identified cavities appears severely discrepant in Figured] 
(note in particular the tightness of the observed correla- 
tion between estimated major and minor axis), which ar- 
gues against a significant population of falsely identified 
cavities. Third, as mentioned above, within the lowest- 
S/N mock subsample (featuring 5,000 net counts), we 
identified cavities in eight out of 18 cases. As can be 
gleaned from Table [2j the corresponding fraction for all 
the real data sets with < 5,000 counts is comparable but 
slightly smaller, seven out of 19. This further supports 
the notion that the incidence of falsely identified cavities 
in our Chandra sample must be small. 

5.3. Surface Brightness Modeling versus Unsharp 
Masking 

The two methods considered here for revealing small- 
scale structure in extended emission each have their 
strengths and weaknesses. The model-fitting approach 
clearly requires data of fairly high signal-to-noise ratio 
to ensure convergence in the fitting process. In addi- 
tion, it assumes a model for the underlying large-scale 
emission which may not always provide a good descrip- 
tion. With unsharp masking, no assumptions are needed 
for the large-scale morphology, but the method requires 
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Fig. 14. — Influence of total source counts, /3, and cavity strength 
on cavity recovery in the residual images from /3-model fitting. 
Top panel shows the detection result for the 27 mock simulations 
with close cavities (D = 15 pixels from the group center) and 
bottom panel the corresponding results for the ones with distant 
cavities (D = 35 pixels). Dark circles represent certain cavity 
detections, gray circles possible ones, and empty circles data sets 
with no detectable cavities. 



a highly non- trivial choice of smoothing scales. Further- 
more, the method normally assumes a circular smoothing 
kernel, so the smoothing on large scales could potentially 
generate spurious cavities in systems with a highly ellip- 
tical X-ray morphology. 

A general conclusion based on all 54 mock images is 
that the model-fitting method seems superior to unsharp 
masking in recovering clearly identifiable cavity struc- 
tures, cf. Figure [13] This is perhaps not surprising, given 
that the underlying group emission in these data is per- 
fectly described by a /3-model (modulo Poisson uncer- 
tainties on source and background counts). However, 
this finding is not limited to our mock data; despite care- 
ful consideration of the smoothing scales employed, Fig- 
ures 2] and Figure [5] show that unsharp masking often 
produces a lower contrast between the cavities and their 
surroundings than the output from the model fitting ap- 
proach. The top panel of Figure [TT] further illustrates 
this point for the mock data, showing another example 
in which cavities in the residual image are more promi- 
nent than those in the quotient image. Another interest- 
ing conclusion can be reached from the bottom panel of 
Figure [HJ which emulates a low-S/N data set with only 
weak cavities. In this case, the two input cavities are 
reasonably well recovered in the residual image, whereas 
the quotient image from unsharp masking would lead 
one to conclude, incorrectly, that cavities are present to 
the east and west of the center of the group emission. 
Hence, even in an idealized case such as this, unsharp 



masking applied to low-S/N data can produce mislead- 
ing results. This serves as a cautionary note if using this 
method alone to search for cavities. Consistency tests 
using other methods, and/or the inclusion of radio data 
to support the existence of cavities, may be crucial for 
reliable identification of cavities in X-ray faint systems. 

In addition, the model fitting approach is occasionally 
able to recover structure that is not seen in the quotient 
images. For example, in the residual image of A1991 
(Figured]), two cavities are clearly visible, while unsharp 
masking only hints at the presence of the northern cav- 
ity. Nevertheless, unsharp masking can still provide a 
very useful supplement to visual inspection and surface 
brightness fitting. This is particularly true in the cases 
where model fitting clearly fails to converge on sensible 
parameters (typically in relatively low-S/N data), and 
for groups in which there are more complicated struc- 
tures which cannot be described by an elliptical /3-model. 
As seen in Figures [4] and [5] a number of our groups belong 
to the latter category. For example, IC 1262 presents a 
sharp surface brightness discontinuity i n the central re- 
gions , possibly related to a cold front (jTrinchieri et al.l 
l2QQ7l ). NGC 4104 shows a clear depression in central sur- 
face brightness, possibly due to cavities being oriented 
along the line of sight. For RXJ 1159+5531, the fitting 
approach fails to clearly uncover the two cavities hinted 
at north-east and south-west of the group center in the 
unprocessed image. These are located at either end of 
an X-ray jet, itself situated at the center of the diffuse 
emission and visible in the raw data (although not clearly 
so in Figure [5}. In the residual image, the cavities are 
overshadowed by some structure generated in the fitting 
process as a consequence of the complex core morphol- 
ogy in this system. In all these cases, unsharp masking 
provides valuable aid in cavity identification. 

6. DISCUSSION 

6.1. Cavity and Central Radio A GN Fractions 

In total, we find that 26 out of our 51 groups (51 ±12%) 
show at least tentative evidence for the presence of cav- 
ities. While this is statistically inconsistent at the ~ 4<j 
level with all our groups showing cavities, our experimen- 
tation with mock images suggests that the real fraction 
might be higher. In particular, since the total number 
of source counts plays such an important role for the de- 
tectability of cavities in our mock data, it is relevant to 
examine the distribution of total source counts for our 
real sample. Figure [15] shows the connection between 
iVtot fo r the observed groups, as defined in Section HJ and 
our cavity detection results. Groups with iVtot > 10,000 
are very likely to be identified as harboring cavities, while 
more than two-thirds of the systems with lower TVtot show 
no clear evidence of these. Successful cavity detection is 
thus clearly biased towards groups with high-S/N data. 
Despite this, there are real (and mock) groups with a 
large number of counts (iVtot > 30,000) for which cavi- 
ties are not detected. In real data, this could certainly 
be because cavities are absent or present a low intrinsic 
X-ray contrast, as suggested by our mock results. Alter- 
natively, it may be due to orientation effects, the impact 
of which we have not tested with the mock data. Also, as 
discussed above, some of our groups show features which 
are clearly not well described by a /3-model, a fact which 
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Fig. 15. — Cavity detectability and the distribution of total 
source counts TV-tot for all 51 observed groups. 



could also complicate cavity recovery when relying on 
surface brightness modeling in the process. 

Nevertheless, given the lower mean S/N of the mock 
systems with non-detected cavities, our inability to de- 
tect cavities in many of the low-S/N Chandra data sets 
may still be consistent with the presence of cavities in 
these systems. As mentioned, for the lowest-S/N mock 
subsample, cavities are identified in 8/18 systems, consis- 
tent with the 7/19 found in real data. Including brighter 
systems (< 15,000 net counts), corresponding values are 
19/36 = 53 ± 15% (mock) and 13/35 = 37 ± 12% (real), 
just consistent at l-cr. The total rate of systems showing 
evidence of identifiable cavities in the mock data of 42/54 
(78 ± 16%) is higher than that of the Chandra sample 
however, but the fractions are still consistent at the 90% 
confidence level. In that sense, our results are broadly 
consistent with the possibility that most, if not all, of our 
real groups do harbor cavities. However, we strongly cau- 
tion against over-interpreting these numbers. Our sim- 
plified mock setup was not designed to reproduce the full 
complexity of the real sample, and the above compari- 
son should be considered indicative at best. Hence, we 
do not claim our results to suggest that all our groups 
contain cavities, only that the actual fraction is likely 
considerably higher than the ~ 50% suggested by our 
analysis. 

A somewhat related question is that of the incidence 
of central radio AGN within our sample. Since six of 
our groups are excluded from such considerations due to 
source confusion or lack of 1.4-GHz coverage, we cannot 
make statistically robust statements for the full sample. 
However, among the remaining 45 systems, 34 (76 ±17%) 
host a detectable central radio source. This fraction is 
higher than the corresponding incidence of detectable 
cavities (23/45, 51 ±13%), but only marginally so, and it 
rises to 21/23 (consistent with 100%) if only considering 
the groups in the combined C- and P-sample. As such, 
the fraction of detectable central radio AGN is compara- 
ble to that of identifiable cavities, especially considering 
that the latter is likely a lower limit to the true frac- 
tion. We do note that while the N-sample is statistically 
underrepresented among groups with detectable radio 
AGN, this is not necessarily the case if adopting a fixed 
cut in Li.4GHz so as to focus on radio-loud AGN alone. 
This is already hinted at in Figure [9j which indicates no 



connection between the presence of detectable cavities 
and the current radio power of the central group galaxy. 
Furth ermore, given that NVSS sensitivity is nearly uni- 
form (jCondon et al.|[l998h and that the weakest detected 
radio source within our sample is also the nearest and has 
a flux just below the NVSS 99% completeness limit, it is 
clear that all our groups could potentially harbor a cen- 
tral AGN with a 1.4-GHz luminosity brighter than this 
source (Li. 4 gh z ~ 1 x 10 20 W Hz" 1 ). 

6.2. Comparison to Results for Clusters 

With half of our sample showing some evidence of 
cavities, our groups are lying between r esults for in- 
dividual giant elliptical galaxies (~ 1/4; iNulsen et al.l 
2009) and the X -ray brightest cool-core clusters (~ 2/3; 
Dunn fc F abian 2006). Our results in Section [5] suggests 
that these differences can be at least partly explained by 
observational biases that favor cavity detection in high- 
flux systems. In addition, as discussed in Section HI the 
surface brightness distribution of groups is generally flat- 
ter than in clusters across the relevant radial ranges. The 
dependence on f3 of the cavity recovery rate (cf. Fig- 
ure [T4j) may therefore also play a role for the above re- 
sults. 

Alternatively, the observed differences in detectable 
cavity fraction may reflect a more fundamental connec- 
tion between the presence of high-contrast X-ray cavities 
and global system properties such as X-ray luminosity or 
total mass. Under the plausible assumption that X-ray 
cavities are generated by a radio outburst in the central 
AGN, the lower cavity fraction in smaller systems could 
be explained by a systematic variation in central AGN 
duty cycle with system mass. If X-ray fainter systems 
generally have shorter central radio outbursts, clearly de- 
tectable cavities are perhaps inflated by the AGN in a 
smaller fraction of those systems. The observed lack of a 
link between the 1.4 GHz luminosity of any radio source 
within the central galaxy and the presence of detectable 
cavities in our sample seems consistent with this possi- 
bility. If the outburst phase is relatively short, bright 
1.4 GHz emission within the central galaxy may already 
have subsided when detectable cavities emerge in the in- 
tragroup medium. This scenario should be testable with 
lower-frequency radio data which trace the older, less 
energetic electron populations from the outburst. Suc h 
studies are currently underway (|Giacintucci et al. 2009). 

With regard to the properties of the detected cavities 
themselves, another finding of our study is that cavity 
radii in the radial and tangential directions are highly 
correlated, with the ratio of the two remaining roughly 
constant regardless of cavity size. The implication of 
equations (P)-© is that bubbles retain a mean axis ra- 
tio of a/b ~ 1.7 as they rise and expand in the sur- 
rounding medium, until eventually shredded by hydro- 
dynamic instabilities. This finding is in excellent quanti- 
tative agreem ent with results for massive galaxy clusters 
(Birz an et al.112004} ). as illustrated in Figure suggesting 
that it is a generic feature of cavity evolution regardless 
of the depth of the local gravitational potential and the 
entropy gradient of the ambient medium (which set the 
characteristic buoyancy rise time of the bubbles). More- 
over, the tight correlation between the size of the cavity 
and the projected cavity distance from the group cen- 
ter (as shown in Figure [8]) indicates that the apparent 
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angle opened by the cavity to the center of the group 
remains nearly constant at ~ 60° during the evolu- 
tion of the cavity. This is also in remarkable agreement 
with the case of clusters (Bir zan et al.ll2004l iDiehl et al.l 
2008), indicating that similar physical processes are re- 
sponsible for bubble evolution, and their disruption, in 
systems covering a large range in total mass. 

6.3. Implications for the Formation and Evolution of 

Cavities 

It is interesting to not e that the gas temperature 
profiles for ou r sample (R asmussen fc Ponmanl |2007; 
iSun et aD2009l ) indicate that 43 of our groups host a cool 
core, with the eight non-cool core systems all belonging 
to the N-sample. Hence, the (C- and P-sample) cavity 
fraction is 60 ± 15% for cool-core groups, in contrast to 
the complete absence of identifiable cavities among our 
non-cool core systems. However, this need not reflect 
any intrinsic differences between cool-core and non-cool 
core systems, since cavity detect ability is conceivably bi- 
ased toward cool-core systems with bright central regions 
(as suggested by Figure [15] and by the increase in cavity 
recovery rate with both iV tot and j3 for our mock data 
sets). On the other hand, barring differences in normal- 
ization, the surface brightness profiles of our different 
subsamples show very similar behavior, with compara- 
ble mean values of /? and r c . As such, there is no clear 
evidence that the surface brightness distribution is more 
centrally peaked for the C- and P-samples, making it 
worthwhile to briefly discuss whether the different cavity 
detection rates for cool-core and non-cool core systems 
reflect a real effect. 

While a detailed discussion of the cool core/non-cool 
core distinction and its possible origins is well beyond 
the scope of this work, we do note that an increased frac- 
tion of detectable cavities among cool-core groups would 
fit the picture of a self-sustained cooling-driven feedback 
loop: Strong central cooling of ICM material feeds the su- 
permassive black hole in the central group galaxy, which 
responds with a radio outburst that creates cavities in 
the surrounding medium and temporarily re- heats it, 
until radiative cooling can re-ignite the process. Some 
support for this comes from the work of IDiehl fc Statlerl 
(2008b), based on a Chandra study of 36 ellipticals, ten 
of which are central group galaxies in our sample (and 
of which eight show evidence of cavities in the present 
work). These authors found a relation between the in- 
ner radial gas temperature gradient and the central radio 
AGN power, in the sense that positive (and steeper) inner 
gradients are preferentially found in optically brighter el- 
lipticals with stronge r ce ntral radio sources (but see also 
Fuka zawa et a l. 2006 and Humphrey et al.ll2006f ). In ad- 
dition, practically all our groups with identifiable cavities 
host a detectable central radio source. Although these 
results suggest a connection between strong central ICM 
cooling, a powerful AGN radio output, and the presence 
of detectable cavities, we emphasize that there is still 
very little direct evidence for central gas cooling to low 
temperatures in groups or clusters. In particular, deep 
X-ray grating spectroscopy of A262, A3581, and HCG 62, 
all belonging to our C-sample, has failed to reveal mate- 
rial cooling to below T ^0.5 keV in their central regions 
([Sanders et al.l 120101 ) . It is also unclear how the cen- 
tral black hole would accomplish the distributed heating 



required to completely halt catastrophic cooling at the 
group core, even if only temporarily. Thus, while the 
above picture may seem appealing, many of its details 
still await observational and theoretical verification. 

Once formed, the properties of cavities themselves 
may help to shed light on the processes responsi- 
ble for their production. In recent years several 
authors have employed numerical techniques to 
study the evolution of buoyant bubbles in galaxy 
groups and clusters, including purel y inviscid hy- 
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Results generally show that although under certain 
circumstances the bubbles could be long-lived, purely 
hydrodynamical evolution usually cannot prevent bub- 
bles from being rapidly shredded by Rayleigh-Taylor, 
Kelvi n-Helmholtz, or Richtmyer-Meshkov instabilities 
(e.g., IDiehl et al.l 120081 ). To prevent bubble disruption, 
additional physics such as magnetic fields, thermal 
conductivity, and viscosity are probably required. With 
certain geome tries, a magnetic fie l d cou l d help to protect 
the bubble (Briiggen & Kaiser 2001; Robinson et al.l 
2004: Uones fc De You ng 2005; Ruszkowski et all 120071: 
Liu et al. 2008; Durs i fc Pfrommerl 120081 : IQ'Neill et"aLl 



2009; Dong & Stone 2009). 



Our observation of cavities at projected distances of 
D > 20 kpc from the group core implies that bubbles are 
reasonably resilient to disruption processes, thus confirm- 
ing the necessity of i ncorporating such additional physics 
in numerical work ([Jones fc De Yo ung 2005). Alterna- 
tively, it has been demonstrated that if a bubble is con- 
tinuously being inflated b y a purely hydrodynamical jet 
(|Pizzolato fc Soker 2001) or a mag netical ly do minated 
jet ([Li et al.l 120061 : iNakamura et al.l I2006L I2007D . under 
some circumstances the development of hydrodynamical 
instabilities could also be suppressed, thus keeping bub- 
bles intact for longer ([Diehl fc Statlerl 120071 ). However, 
the fact that we find no clear link between the presence 
of detectable cavities and the current luminosity of any 
central radio source seems to argue against this mecha- 
nism as a general feature in our groups. 

To briefly explore the possibility of distinguishing be- 
tween different model predic tions for cavity ev olution, 
we can adopt the approach of IDiehl et al.l (|2008f ). These 
authors consider a range of different models for the for- 
mation and evolution of cavities, under the reasonable as- 
sumptions that the bubbles rise subsonically in the ICM 
in pressure equilibrium with their surroundings, with 
the ICM background density described by a /3-model. 
The models considered include one of purely hydrody- 
namic, adiabatically expanding bubbles containing gas 
with an adiabatic in dex T = 5/3 (here denoted "AD53" 
for consistency with IDiehl et al.l 12008). and a model in 
which bubbles are generated by a current-dominated jet 
("CDJ", see, e.g.. iLi et al.ll2006f ) and have magnetically 
dominated pressure. In both these models, predicted 
bubble radii a grow with clustercentric distance D as 
a oc [1 + (D/r c ) 2 ] a , with a = P/(2T) and a = 3/3/4 for 
the AD53 and CDJ models, respectively. The asymp- 
totic growth of the bubbles (at D ^> r c ) is thus a pow- 
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Fig. 16. — Cavity major axis a as a function of projected group- 
centric distance D normalized by r c . Dotted line ("AD53") shows 
the predicted evolution of a purely hydrodynamical bubble, and 
dashed line ("CDJ") that of a magnetically dominated bubble in- 
flated by a current-dominated jet (see text for details), both for an 
arbitrarily chosen initial bubble size of a = 3 kpc. 

erful discriminator of these diffe r ent m odels. Among the 
models discussed by lDiehl et al.l ((2008), AD53 represents 
the slowest and CDJ the fast est expansion of th e bubbles. 
When plotted against £>/rv. lDiehl et al.l (|2008l . their fig- 
ure 2) find that observed cavity sizes at D > 4r c suggest 
that bubbles expand relatively fast, with CDJ the pre- 
ferred model. 

For comparison, we show a similar plot for our cav- 
ity sample in Figure HU Overplotted are model predic- 
tions for the AD53 and CDJ models, both assuming a 
fixed value of j3 = 0.45, the mean value for our groups 
with detectable cavities. It is clear that we cannot eas- 
ily distinguish between these two simple models. This 
is largely due to a dearth of system s with dista n t cav - 
ities in our sample, since — unlike iDiehl et al.l ([2008) 
— we have only a handful of cavities detected well out- 
side the ICM core radius. A larger group sample, par- 
ticularly including systems with distant cavities, would 
be required to make robust statements on the nature of 
cavities in groups. However, since the cavities within our 
sample appear to obey the same scaling relations as seen 
in more massive systems, constrain ts from studies of cav- 
ities in clusters (jDiehl et al.l 12008) are likely to apply in 
our groups as well. 

7. CONCLUSIONS 

Based on Chandra archival data of 51 galaxy groups, 
we have presented the first systematic search for, and 
analysis of, X-ray cavities in a large sample of groups. 
At least nine of these groups have previously been re- 
ported to show evidence of such cavities. In the present 
work, the cavities are identified from Chandra data alone, 
either by subtracting an elliptical /3-model fit from the 
X-ray surface brightness distribution, or by performing 
unsharp masking using two Gaussian smoothing kernels. 
Based on visual inspection of the resulting images, 13 
groups (~ 1/4) within our sample are identified as clearly 
harboring cavities, and another 13 groups (~ 1/4) show 



tentative evidence for cavities. All of these 26 systems 
show clear evidence for a cool core. The remaining 25 
groups have a lower cool-core fraction (17/25 ~ 70%) 
and do not show convincing evidence for the presence 
of cavities. The total fraction of groups with some evi- 
dence of cavities in our sample is thus 50%, rising to 60% 
if only considering the cool-core systems. In this respect, 
our groups are lying betw een results for indiv idual giant 
elliptical galaxies ( ~ 1/4; Nulsen et al. 2009|) and cool- 
core clusters (~ 2 /3: lDunn fc F abian 2006). We found no 
clear link between the current 1.4 GHz radio luminosity 
of the central brightest group galaxy and the detection 
of cavities. 

To test our ability to identify cavities using the two 
adopted approaches, we generated a set of mock im- 
ages designed to mimic typical Chandra images of the 
observed groups. For these data, we find that sur- 
face brightness modeling is generally superior to unsharp 
masking in reliably recovering cavity properties, a con- 
clusion echoing one already hinted at from our analysis of 
the real data. The results also show that cavity recovery 
is strongly influenced by a number of factors, including 
the projected distance between the cavity center and the 
group center, the significance of the surface brightness 
depression represented by the cavity, and most impor- 
tantly, the signal-to-noise ratio of the diffuse group emis- 
sion. Cavity detectability is highly biased toward X-ray 
bright systems, implying that the observed cavity frac- 
tion in our sample is likely a lower limit. 

We find tight correlations between the radial and tan- 
gential radii of the cavities in our group sample, and 
between the size and projected groupcentric distance of 
cavities. While we are not able to clearly distinguish be- 
tween simple model predictions for the evolution of cav- 
ities on the basis of our sample alone, we note that the 
observed correlations are in excellent quantitative agree- 
ment with results for more massive clusters. This sug- 
gests that very similar physical processes are responsible 
for cavity evolution and disruption in systems covering a 
large range in mass. 
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a Chandra Observation ID. 



Cleaned Chandra exposure time. 
c Luminosity distance from the NASA/IPAC Extragalactic Database (NED), in 
the reference frame defined by the cosmic microwave background. 
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TABLE 2 
Fit results and cavity properties 
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TABLE 2 
Fit results and cavity properties 

UGC842 N 3592 22.23 

UGC2755 N 0.679±J5;?f? 408 24.24 

a Classification according to whether the group contains certain cavities (C), possible cavities (P), or no detectable cavities (N). 

k Projected distance from cavity center to the X-ray center of the group. 

c Cavity size in tangential direction (projected major axis a defined in Section I3.3J) . 

^ Cavity size in radial direction (projected minor axis b defined in Section |3.3|) . 

6 Number of 0.3—2 keV photons from diffuse emission on the central CCD. 

Radio luminosity at 1.4 GHz of any radio source within the central brightest group galaxy, extracted from the NRAO VLA Sky Survey (NVSS; 
iCondon et al.lfl998T) . See text for details. 

^ The cavity in NGC 4104 is at the very center of the X-ray emission. 



